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^ ABSTRACT 

Galaxy clusters grow by mergers with other clusters and galaxy groups. These mergers 
r— i create shock waves within the intracluster medium (ICM) that can accelerate particles 

to extreme energies. In the presence of magnetic fields, relativistic electrons form 
large regions emitting synchrotron radiation, so-called radio relics. Behind the shock 
front, synchrotron and inverse Compton (IC) losses cause the radio spectral index to 
steepen away from the shock front. In a binary cluster merger, two shock waves are 
generated which move diametrically outwards along the merger axis. Two radio relics 
Q can then form on both sides of the cluster center. An example of such a cluster is 

CIZA J2242. 8+5301, where very clear spectral steepening in the downstream region is 
c/) observed. The main relic has a total extent of 1700 kpc, while its width is only 55 kpc. 

Together with the high observed polarization fraction, this implies the relic is seen 
very close to edge-on which makes it easier to constrain the merger geometry. Here 
we present hydrodynamical simulations of idealized binary cluster mergers with the 
aim of constraining the merger scenario for this cluster. From our simulations, we find 
OO that CIZA J2242. 8+5301 is probably undergoing a merger in the plane of the sky (less 

0^ then 10° from edge-on) with a mass ratio (Mi : M?) of about 2:1, and an impact 

C*~) parameter < 400 kpc. We find that the core passage of the clusters happened about 

1 Gyr ago. We conclude that double relics relics can set constraints on the mass ratios, 
QO impact parameters, timescales, and viewing geometry of binary cluster mergers, which 

C*^ is particularly useful when detailed X-ray observations are not available. In addition, 

t-H the presence of large radio relics can be used to constrain the degree of clumping in 

the outskirts of the ICM, which is important to constrain the baryon fraction, density 
and entropy profiles, around the virial radius and beyond. We find that the amplitude 
• i-H of density fluctuations, with sizes of < 200 kpc, in the relic in CIZA J2242. 8+5301 is 

^ not larger than 30%. 

H 
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1 INTRODUCTION 

Radio relics are diffuse, steep-spectrum radio sources found 
in merging galaxy clusters. They have been divided into 
three groups (Kempner et al. 2004 [): radio gischt, radio 



phoenix and AGN relics. Here we will only be concerned 
with so-called radio gischt relics, these are elongated arc-like 
radio sources with sizes of up to 2 Mpc (e.g., |Ensslin et aTj 
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1998 1. They are mostly found in the outskirts of galaxy clus- 
ters. Recent observations have given support to the hypoth- 
esis that they trace shock fronts in which particles are accel- 
erated via the diffusive shock acceleration mechanism (DSA 
Krymskii|[l977l |Axford et al.||1977| |Bell||1978a |b| |Blandford| 
fc Ostrikerl |1978[ |Drury||1983| |Blandford fc' 



Jones fc Eliison||1991| |Malkov fc O'C Drury|2001 1 



Eichler][1987 



rd 

57; 



A particularly interesting "class" of radio gischt are 
the so-called double-relics. In this case two relics are dia- 
metrically located on both sides of the cluster center (e.g., 
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Bonafede et al.||2009| |van Weeren et al.||2009| |Venturi et al. 


2007| |Bagchi et al.|2006; Rottgering et al.|1997 


van Weeren 


et al.|2010||Brown et al.|2011 van Weeren et al. 


20TT1 Bagchi 


et al. 2011|). These double relics are very rare, with only 



about ten known so far. However, most of them were found 
in the last few years suggesting that more of them await 
discovery, especially with future deep radio surveys. Large 
cosmological simulations that include radio emission from 
shocks, also indicate that double radio relics are more com- 
mon (e.g., |Skillman et al.|2011 Battaglia et al.||2009 Hoeft 



et al.|2008| |." 

Double radio relics are thought to trace outward mov- 
ing shock waves from a binary cluster merger, which develop 
after core passage of the two subclusters. In the DSA sce- 
nario, synchrotron and IC losses cool the particles at the 
back of the shock front. If the relics are seen roughly edge- 
on this creates a radio spectral index gradient towards the 
cluster center. The line connecting the two relics should rep- 
resent the (projected) merger axis of the system. Thus, dou- 
ble relics in principle can be used to set constraints on the 
merger timescale, mass ratio, impact parameter and viewing 
angle. 

An example of a cluster hosting a double radio relic 
is CIZA J2242. 8+5301 ( |van Weeren et al.||2010| ). The clus- 
ter is located at 2 = 0.1921 and has an X-ray l uminosity 



of 6.8 x 10 44 ergs -1 , between 0.1 and 2.4 kcV (Kocevski 



et al. 20071. With this X-ray lumino sity the total mas s of 
the cluster is roughly 5 .5 x 10 14 M Q ( |Pratt et aL|2009| ). In 
van Weeren et al. (2010 1 we presented Westerbork Synthesis 



Radio Telescope (WSRT) and Metrewave Radio Telescope 
(GMRT) observations of this cluster that showed the pres- 
ence of an impressive double relic system. In the north of 
the cluster, a large radio relic is found with a total extent 
of 1700 kpc (we did not include the patch of emission at 
RA 22 h 42 m 20 s , Dec +53° 08' 45"). To the south of the 
cluster center a second fainter relic is found with an extent 
of 1450 kpc. The relics are located along the major axis of 
the elongated X-ray emission from the ICM, while the relics 
themselves are orientated perpendicular to this axis. This 
configuration is expected for a "clean" (i.e., a merger event 
without much substructure) binary merger event in which 
shock waves are propagating diametrically outwards to the 
north and south. Within the shock particles are proposed to 
be accelerated by DSA. The relics are separated by a dis- 
tance of 2.8 Mpc. The northern relic is strongly polarized 
at 4.9 GHz, with a polarisation fraction of ~ 50%. The po- 
larization magnetic field vectors (corrected for the effect of 
Faraday rotation) are parallel to the major axis of the relic. 
The magnetic fields are thus mainly orientated within the 
plane of the shock wave and the relic must be seen close to 
edge-on ( |Ensslin et al.|1998[ |. We also observe clear steepen- 
ing of the radio spectrum (with the spectral index changing 
from —0.6 ~ —0.8 to more than —2.0 between 0.61 and 
2.3 GHz) in the direction of the cluster center, across the 
width of the relic. Between the two relics additional diffuse 
emission is found, which is also seen in a deep 150 MHz 
GMRT image ( |van Weeren et al.|20f7 i. 

However, there are still quite a few puzzles that sur- 
round radio relics, such as: 

• Which processes accelerate electrons so efficiently at rel- 
atively low Mach number (M ~ 2 — 4) shocks? 



• What produces the magnetic fields inside relics? 

• Why do some relics have very sharp edges while others 
appear very fuzzy? 

• Under which conditions do relics form? When do we see 
single and when double relics? 

Many of these questions touch physical processes that 
are poorly understood, such as diffusive shock acceleration 
and magnetogenesis in collisionless shocks, and relate to re- 
gions of the cosmos of which we know very little. At the same 
time, outskirts of galaxy clusters have attracted consider- 
able interest for cosmology as they are governed primarily 
by simple gas physics and gravity and much less by com- 
plicated processes such as radiative cooling, star formation 
and AGN feedback. However, information about the periph- 
ery of galaxy clusters is scarce. Recent X-ray observations 
with Suzaku have measured surface brightness profiles out 



to the virial radius for a number of nearby clusters ( 


Bautz 


et al.|2009||George et al.|2009| 


Rciprich et al.|2009||Hoshino 


et al.] 


2010| [Simionescu et al. 


2011 1. The results indicate 



significant deviations from extrapolated profiles of the gas 
mass fraction which are explained by an increased dumpi- 
ness of the ICM that occurs on scales of several 100 kpc or 
smaller. As the cluster outskirts are difficult to access via X- 
ray observations owing to the low surface brightness at large 
cluster-centric radii, important information may be derived 
from radio relics which illuminate the outskirts of merging 
galaxy clusters. 

In this paper we will describe simulations with the 
aim of reconstructing the merger event of the cluster 
CIZA J2242. 8+5301. In particular we will focus on the mass 
ratio and impact parameter of the merger and how they can 
be constrained from the observations and simulations. We 
further attempt to infer properties of the ICM at distances 
beyond the virial radius from the observed features of the 
radio relics, such as the morphology of the emission. 

In Sect. [2] we will shortly discuss previous work on ide- 
alized models of cluster mergers, and in Sect. [3] we describe 
the method and adopted initial conditions. This is followed 
the results in Sects. [4] and [5] We end with a discussion and 
summary m Sects. [6] and [7] 



2 OVERVIEW: SIMULATIONS OF GALAXY 
CLUSTER MERGERS 

Numerical simulations of idealized binary galaxy cluster 
mergers were performed by Roetti ger et al.|p993[);|Schindler 
& Mueller|(|1993|l;|Burns et al.l (119941); iPearce et all (119941); 



Roettiger et al.| (|1997[); |Ricker ( |1998[ ); |Roettiger fc Flores 
(2000) JTRlcker~fc Sarazinl ( |200l| 



Roettiger et al. ( 1993[ ) simulated a head-on merger with 



a mass ratio of 8. The clusters were modeled as King-spheres. 
They used a hybrid Hydro/N-body code, simulating the 
hydrodynamical component with the ZEUS-3D code on a 
nonuniform grid, and the N-body component with the Hern- 



quist treecode (Hernquist 19871. Their simulation showed 



the development of a single shock that reached a Mach num- 
ber (M) of more than 4. In addition the core of the main 
cluster was heated to 3 x 10 8 K. The simulated X-ray sur- 
face brightness was double peaked and elongated parallel to 



the merger axis. Burns et al. ( 1994 1 simulated the collision 
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between the Coma cluster and the NGC 4839 galaxy group during cluster mergers and Akahori & Yoshikawa (20101 in- 



using the same code as described by Roettiger et al. ( 1993 1 



Schindler & Mueller ( 1993 I simulated collisions of clus- 



ters using a uniform grid-PPM (piecewise- parabolic method) 
code for the gas dynamics, while an N-body code was used 
for the collisionless component consisting of dark matter and 
galaxies. They showed the development of outward propa- 
gating shock waves generated during the collisions between 
subclusters. 

None of the works mentioned above included the bary- 
onic contribution to the gravitational potential. Head-on 
equal mass collisions between clusters were simulated by 
|Pearce et al.| ( |1994[ ) using smoothed-particle hydrodynam- 



ics (SPH), including both baryons and dark matter. Roet 



tiger et al. (19971 used a PPM/particle-mesh code to simu- 



late a collision intended to resemble the cluster Abell 754. 



Roettiger et al. ( 1999 1 presented the first three-dimensional 



numerical magnetohydrodynamical (MHD) simulations of 
the magnetic field evolution in merging clusters of galaxies, 
which they subsequently applied to model the double radio 
relics in Abell 3367 ( |Roettiger et al.|1999| |. They concluded 
that the radio structures arise from a slightly off-axis merger 
that occurred nearly in the plane of the sky approximately 
1 Gyr ago, with the subcluster having a total mass of 20% 
of the primary cluster. |Roettiger fc Flor cs ( 2000 ) also mod- 
eled the X-ray substructure seen in Abell 3266 as an off-axis 
merger with a mass ratio of 2.5 : 1. 

|Ricker| ( |1998[ ) simulated off-center collisions of equal 
mass clusters using a PPM code and an isolated multigrid 
potential solver. They simulated off-center collisions with 
0, 5, and 10 times the cluster core radius and studied how 
the virialisation time, X-ray luminosity, and structure of the 
merger depended on the impact parameter of the collision. 
They found that the increase in X-ray luminosity due to the 
merger depends on the impact parameter of the collision. 



Ritchie & Thomas ( 2002 1 also reported that mergers can 



lead to large increases in bolometric X-ray luminosities and 
emission-weighted temperatures. Cool cores are completely 
disrupted by equal mass head-on mergers, while for mergers 
with mass ratios of 8 : 1 the cooling flow restarts within a 
few Gyr. For mergers with impact parameters of 500 kpc 



the cool-core is not disrupted. Ricker & Sarazin ( 2001 1 in- 



vestigated cluster collisions, varying both the mass ratio and 
impact parameters. They found that merger events created 
large-scale turbulent motions with eddy sizes up to several 
hundred kiloparsecs. 

FLASH dFryxell et al.l 120001) high-resolution N- 



body/hydrodynamics simulations have been carried out by 
|Zu Honeet~aLl ( |2009| ) to model the cluster CI 0024+17, a pro- 
posed merger of two clusters, with the interaction occurring 
along our line of sight. Also, gas sloshing initiated by mergers 
with subclusters, and the effects of mergers on the entropy of 



the ICM were investigated by ZuHone et al. (20101; ZuHone 



(2011 1. Various simulations of cluster merger were also pre- 
sented by[Takizawa| |T999| [2000]) ; |Takizawa fc Naito| fMity ; 
Takizawaj ( |2008[ ). In |Takizawa| ( |2008| l N-body + magneto- 



hydrodynamical simulations were performed from which is 
was found that cluster mergers cause various characteris- 
tic magnetic field structures. The magnetic field component 
perpendicular to the collision axis is amplified by the merger 
which results in a cool region wrapped by field lines. |Mc-| 
Carthy et al. ( 2007 1 focussed on the generation of entropy 



vestigated how the ionization equilibrium state is affected 
by cluster mergers. They report that the ICM significantly 
departs from the ionization equilibrium state at the location 
of shocks. 



Poole et al. (2006) carried out a suite of SPH simula- 



tions of merger of clusters whose initial conditions resembled 
relaxed cool core clusters. They investigated the efficacy of 
(i) centroid variance, (ii) power ratios, and (iii) X-ray sur- 
face brightness/projected mass displacement to quantify the 
degree to which observed clusters are disturbed. From this it 
was found that the centroid variance gives the best measure 
of the state of clusters and how far they are from virial and 
hydrostatic equilibrium. 

Cluster mergers can decouple the baryonic matter com- 
ponent from the dark matter (DM) which causes and off- 
set between the gravitational center (measured from lens- 
ing) and X-ray center of the cluster. This was first observed 
for the "Bullet cluster" (1E0657-56, |Clowe et ai^|2006 |. 
|Springel fc Farrar] ( |2007[ ) presented hydrodynamical models 
of galaxy cluster mergers to reproduce the dynamical state 
and mass models (from gravitational lensing) of the "Bul- 
let" cluster (1E0657-56). They showed that the size of the 
spatial offset between the baryonic peak and mass is quite 
sensitive to the structural details of the merging systems. 
This offset was earlier also found in numerical models by 
iTakizawal |2006| l ~ 



Mastropietro & Burkert (120081) presented 



detailed N-body/SPH simulations of the system. They found 
that the X-ray morphology is best simulated with a mass ra- 
tio of 6 : 1 and an impact parameter of 150 kpc. 



3 NUMERICAL METHOD 



We used the FLASH 3.2 framework ( |Fryxell et al.|2000[ ) for 
the simulations which includes standard hydrodynamics and 
gravity. We did not include the effects of radiative cooling, 
as this is not very relevant for the spatial and temporal time 
scales considered here. Furthermore, the hydrogen number 
density was assumed to be related to the electron number 
density as nn = 0.6n o . We chose outflow boundary condi- 
tions. We simulated a box with a size of 5 x 5 x 5 Mpc. A 
maximum of 6 refinement levels are allowed, resulting in a 
maximum resolution of 3.25 kpc. The minimum refinement 
level was set to 3 giving a resolution of 26 kpc. The simu- 
lations were run at this relatively low-resolution in order to 
study a larger range of merger parameters. 

For the gravitational potential of each subcluster we 
assume hydrostatic equilibrium and spherical symmetry, 



dr 



T(r) dp(r) dT{r) 



+ 



p(r) dlogr dlogr 



(1) 



where T(r) and p(r) are the radial density and temperature 
profiles, ttih the hydrogen mass, and p, the mean molecu- 
lar weight (we adopt p = 0.6). We move the center of the 
gravitational potential of the merging subcluster around the 
fixed potential of the main cluster, ignoring the interactions 
between the dark matter. 
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3.1 Radio emission from shocks 

We use passive tracer particles in the simulations to model 
the radio emission. At the start of the simulation the parti- 
cles are distributed proportional to the density throughout 
the complete computation volume. The advantage is that in 
this way each passive tracer particle represent roughly the 
same amount of mass which simplifies the normalization of 
the radio emission, as the contribution from every particle 
can be simply added up. For each tracer particle, 5 million in 
total, we record the position, velocity, density and temper- 
ature as a function of time, writing the results to disk every 
30 x 10 B yr. In order to create synthetic radio maps, we lo- 
cate all particles that have passed through the shock waves 
in the simulation by looking at the variation of the temper- 
ature as function of time for each tracer particle. We define 
a particle to have passed through a shock if the temperature 
between two successive outputs spaced 30 x 10 6 yr apart in- 
creases by more than a factor of 1.5 (no radio emission is 
generated below this number). For each particle that passes 
through the shock we then compute the Mach number (M, 
using Rankine-Hugoniot jump conditions and adiabatic in- 
dex 7 = 5/3 ( |Landau fe Lifshitz|p59] |, compression ratio, 
entropy ratio, and the time since it has passed through the 
shock. With these parameters we compute the synchrotron 
emission at a given frequency v using the JP model ( jJaffe fc| 
Perola|1973 1 , taking into account the spectral ageing due to 



synchrotron and IC losses. We apply eq. 5 and 6 f rom|Komis- 
sarov fc Gubanov] (|1994[> a nd use to = 0, see also IKardashev 
1 1962[); |Pacholczyk| (| 1970 1 . 



The amount of synchrotron emission for each tracer par- 
ticle is recorded. A radio map is then simply computed by 
integrating the radio emission from each tracer particle in 
the computational volume along a chosen line of sight. 

We take Bomb = 4.6 pGauss the equivalent magnetic 
field strength of the CMB at z = 0.1921. To compute the 
radio emission the magnetic field is needed. Since we do not 
include magnetic fields in our simulation we adopt constant 
values for the magnetic field B. The spectra are then nor- 
malized using the method described by |Hoeft fc Briiggen 
(20071, which takes into account the efficiency of acceler- 



ation as function of Mach number. We start with a spec- 
trum of suprathermal electrons with a power-law distribu- 
tion n(E) oc E~ s , with E the energy. The index s, with 
s = 1 — 2ai n j, and Qinj the injection radio spectral index, 
depends only on the compression ratio (or Mach number), 
r, of the shock front: 



r + 2 

S= ^T 

and (e.g., Giacintucci et al.[2 008) 

1 _ M 2 + 1 
Qinj - 2 _ M 2 -l 



(2) 



(3) 



We assume that a small fraction £ c of t he thermal en- 
ergy goes into the acceleration of electrons (Keshet et al. 



|2004[ ). Furthermore, we take into account that the electrons 
can only be acceler ated to a finite maximum energy (-E max ) 
( jKeshet et al. |2003 1 . The electron spectrum with a smoothed 
high-energy cut-off is then given by 



dn 
di? 



n C C s p cc 2 



6 ^max 

elsewhere 



(D 



where e = E/rn e c 2 . The normalization constant C ap ec gives 
the fraction of electrons at E = m e c 2 . 

Following |Hoeft fc Briiggen (20071, we assume that the 
thermal Maxwell-Boltzmann distribution goes over continu- 
ously into the power-law spectrum of the suprathermal elec- 
trons. Hoeft & Briiggen (2007) show that the normalization 



of the electron spectrum is given by 



Ud m p (q — 1) 1 



(5) 



where Ud is the downstream velocity, q the entropy ratio 
(Sd/Su) across the shock, and 



^soec — 



de e 1 



(6) 



The minimum energy, e m i n = _B m i n /m c c 2 , is the energy 
above which electrons are considered to be suprathermal. 
This results is an implicit equation for -E m m that has to be 
solved simultaneously with the normalization of the spec- 
trum. We directly apply the code from |Hoeft fe Briiggen| for 
this. 



3.2 Initial conditions 

For the initial conditions we start with two spherically sym- 
metric (sub)clusters with masses Mi and M2 separated by 
a distance d. Subscript 1 always refers to the more mas- 
sive cluster (i.e., Mi > M2). The density profiles of both 



subclusters are described by single ,/3-modcls (Cavalicrc & 



Fusco-Femiano 1976) 



p(r) 



Po 



1+ ( - 



(7) 



where po the central density, r, 
to fix to 2/3 



10" 26 g cm - 



(e.g. 



Bohringer et al 
3 (0.02 particles per cm 13 ' 



the core radius. We decided 
|2000[ ) and po = 2 x 
for both subclusters. 



This value for po is roughly the average of the values found 
by Croston et al. (2008) for a sample of 31 clusters. They 



found that the central densities did not show any clear trend 
with the overall mass and temperature of the clusters. 

We used the M — T scaling relation from |Mantz et al.| 



(20101 to get the global average temperature from the clus- 



ter s mass. 



1Q 0. 88+0.48 log 1Q (M/1O 14 M ) rj^yi 



(8) 



For the radial temperature distribution we use a polytropic 
profile (e.g., |Pratt fc Arnaud|2002"l l 

T(r) = T avg rp(r)^" 1 (9) 

with 7" — 1.1 and 7t = 1.2 and T avg given by Eq. [8] In ad- 
dition we allow for a temperature drop in the center of the 
cluster to simulate a relaxed cool core cluster (e.g., [Fabian 
et al.||1991||Peterson fc Fabian||2006[ ). T he temperature de- 



cline is the central region is given by (Allen et al. 
Vikhlinin et al.||2006[ ) 



2001 



T cool (r) = 



x + l 



^"cool 



(10) 



where T m i n /T avg is the relative temperature drop in the cen- 
ter, r coo i the radius of the cool core, and a CO oi controls de rate 
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of the temperature decrease. The final temperature model 
is given by 



T fu n(r) = T(r)T cool (r) 



(11) 



For the merger scenario we follow Sarazin ( 2002 1, where 
two subclusters merge at time i me r g e (the age of the Universe 
at the time of the merger). The two clusters have a large 
initial separation, do - The value of do does not strongly affect 
the merger because the infall velocity approaches free-fall 
from infinity. If we assume that the clusters dominate the 
mass in the region of the Universe they occupy 



do w 4.5 



Mj + M 2 
1O 15 M 



1/3 



tn 



10 1 



yr 



2/3 



[Mpc] . (12) 



For CIZA J2242. 8+5301 (z = 0.1921) i merg c is 11.0 Gyr. 

To speed up the simulations the centers of the two sub- 
clusters are placed at a distance d = 2 Mpc from each other. 
This we define as t — 0. The relative vel ocity between the 
clusters at that starting point is given by Sarazin ( 2002 1 as 

\ 1/2 

\ / J X -1/2 

[km s _1 ] 



2930 



Mi - 



M 2 



_d_ 

'hi 



M 



1- 6 



d 

1 Mpc 



(13) 

where b is the impact parameter of the collision. Using the 
starting location, initial velocities and masses of the subclus- 
ters we move the center of the gravitational potential of the 
merging subcluster (M2) around the fixed potential of the 
"main" cluster (Mi), ignoring the interactions between the 
dark matter and treating the clusters as point masses. This 
fixes the subcluster's orbit around the main cluster. The 
point (a;, y, z) = (0, 0, 0) is defined by the center of mass of 
the system. The gravitational acceleration of the gas was 
then computed using Eq. [I] 

To reconstruct the merger that gave rise to the double 
radio relic CIZA J2242. 8+5301, we varied the mass ratio and 
impact parameter in the simulation. In the post-processing 
(to create the radio maps) we also varied the viewing angle. 
The total mass of the two (sub)clusters (Mi +M2) was kept 
constant to 5.5 x 1O 14 M0. The values for r c are fixed by 
the total mass of the system and the required mass ratio 
(these are reported in Table [T]). In Sect. [6] we will discuss 
the effect on the simulations if we take different values for 
/3 and r c . For the default runs we chose T m i n /T aV g = 1, i.e., 
no temperature drop in the cluster center. In Sect. |4.1.5| 
we discuss the effect of including a central cool core to the 
most massive subcluster. For all runs we place the merger 
axis along the y-axis with the least massive cluster located 
at positive y-value and the initial velocity (given by Eq. 13 
being in negative y-direction) . 



4 RESULTS 

The names and parameters for the different FLASH runs 
are given in Table [I] We start with a description of the 
basic hydrodynamical quantities (the density, temperature, 
and the velocity field) from the simulations. We will discuss 
the results for a merger with a mass ratio of 2 : 1 and zero 
impact parameter (R21b0; see Fig. [TJ , and for a merger with 
the same mass ratio but with an impact parameter of 4 
(4r c ,i; R21b4), see Fig. [2] 



The R21b0 merger starts with the subcluster approach- 
ing the main cluster. Two shocks, with M ~ 1.15 at 
t — 600 Myr, develop as the cores approach each other. After 
core passage, at t = 1.25 Gyr, a pair of stronger symmetrical 
shock waves form that travel diametrically outwards in front 
of the two subclusters. The shock waves increase in size as 
the time advances, although their Mach numbers decrease. 
At t = 1.95 Gyr the maximum Mach numbers are about 3.1 
and 2.3 for the top and bottom shock waves, respectively. 
At t = 2.28 Gyr the two shock are separated by 2.8 Mpc, as 
is the case for the double relics in CIZA J2242. 8+5301. The 
Mach numbers at this time have decreased to about 2.7 and 
2.2, for the top and bottom relics, respectively. Thus, the 
most massive cluster with the largest core radius develops 
the largest (and slightly stronger) shock wave in front of it 
after core passage. The velocity field at t — 2.3 Gyr is shown 
in Fig. [3] (left panel). 

Mergers with different mass ratios, ranging between 1 
to 5, show a similar picture. The size and strength of the 
shock waves depend on the mass ratio (which in turn deter- 
mines the core radii of the subclusters). In collisions with a 
larger mass ratio the shock waves in front of the less massive 
subcluster reduces in relative size and strength. 

A 2 : 1 merger with an impact parameter of b = 4r Ci i 
(R21b4) develops two asymmetric spiral-like shock waves, 
and again the two shocks have different sizes depending 
on the mass ratio used. The amount of asymmetry deep- 
ens on the impact parameter. The Mach numbers across 
these asymmetric shocks decrease more slowly on the left 
side of the merger axis. At t = 2.1 Gyr the maximum Mach 
number is about 2.9 and 2.2 for the top and bottom shock 
waves, respectively. This decreases to about 2.7 and 2.1 at 
t — 2.4 Gyr when the two shocks are separated by 2.8 Mpc. 
The velocity field at t — 2.4 Gyr is shown in Fig. [3] (right 
panel) . 

4.1 Radio maps 

In this section we present the simulated radio maps for a 
range of mass ratios, impact parameters, viewing angles 
and radial density/temperature profiles, and compare them 
to the observed WSRT radio map at 1382 MHz. To create 
the radio maps, we integrated the radio emission from each 
tracer particle that passed through the shock along the line 
of sight. All radio maps are normalized to the peak flux in 
the simulated map and are smoothed with a 2-dimensional 
Gaussian of 50 x 40 kpc to match the resolution of the WSRT 
image. This smoothing is also needed as the radio emission 
is computed using a discrete number of particles with in- 
sufficient sampling to create images at the resolution of the 
hydro-grid. For making the radio maps, we assumed a con- 
stant magnetic field strength of B — 5.0 ^iGauss, see [van] 



Weeren et al. (2010) 



It takes about 1.0 Gyr for the two relics to be separated 
by 2.8 Mpc after core passage. This number varies very little 
for mass ratios between 1.5:1 and 3:1. For impact parame- 
ters b > kpc the time since cores passage only increases 
slightly to about 1.1 Gr for b = 5r c .i (e.g., Fig. [3|. There- 
fore, the time since core passage does not depend much on 
the adopted initial conditions such as the mass ratio and im- 
pact parameter. However, the time since core passage will 
also depend on the adopted total mass of the system and on 
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Table 2. Relic sizes 





top relic 


bottom relic 


CIZA J2242.8+5301 


1700 


1450 


run R1.51b0 


1730 


1570 


run R21b0 


1670 


1220 


run R2.51b0 


1680 


930 


run R31b0 


1690 


660 



the temperature profile, since this affects the sound speed. 
Detailed X-ray observations would be helpful here to bet- 
ter constrain the total mass and measure the temperature 
profile. 

4-1.1 Mass ratio 

We varied the mass ratios for the collisions between 1 : 1 
and 3:1. The mass ratio (which also set the core radii) 
affects the size of the two outwards moving shock waves after 
core passage. The resulting radio maps for mergers with four 
different mass ratios are displayed in Fig. [4] and the largest 
linear extent of the relics are listed in Table [2] The radio 
maps display a double relic system, with the front of the two 
relics tracing the shock waves. Towards the cluster center 
and the left/right edges of the relics, the surface brightness 
decreases. The size of the shock waves (the A4 > 2 part) 
directly reflects the extent of the relics in the radio maps. 

We compare the synthetic radio maps to the observed 
WSRT map at 1382 MHz. The observed and synthetic 
maps are aligned by eye to the brightest northern relic in 
CIZA J2242. 8+5301, by simple translation and rotation. 
With a visual comparison and the taking into account the 
sizes listed in Table. [2j we find the best match between the 
observed and simulated radio maps for a mass ratio around 
2 : 1 to 1.5 : 1 For more equal mass ratios the bottom relic 
increases to a size that is larger than the observed size. The 
opposite happens for larger mass ratios. When we compare 
the simulated radio images with the temperature map in 
Fig. [3] (left panel) we note that the relic emission does not 
trace the full extent of the shock waves. This is a conse- 
quence of the non-linear relation between the Mach num- 
ber and normalization of the radio spectrum, discussed in 
Sect. 13. ll Below Ml ~ 2 there is no observable radio emission 
generated. Another reason for the gradual disappearance of 
the radio emission at the right/left edges of the relics is that 
the amount of projected shock surface along the line of sight 
decreases. 



4-1.2 Impact parameter 

Taking the 2 : 1 merger event, we varied the impact param- 
eter, see Fig. [4] As discussed in the beginning of Sect. [4j 
two asymmetric shock waves develop, breaking the cylin- 
der symmetry (or reflection symmetry in the resulting im- 
age) around the merger axis. The variation in Mach number 
across the shock front and the amount of projected shock 
surface along the line of sight cause the radio emission to 
fade away more slowly on one side of the radio relics. Inte- 
grated brightness profiles for the top relic across the y-axis 



are shown in Fig. [5] (top panel) . The total extent of the relic 
decreases by about 10% from b — kpc to b — 5r c ,i. 

The CIZA J2242. 8+5301 relics are more or less symmet- 
rically located along the proposed merger axis, although the 
brightness profile at the east end of the northern relic is not 
well determined because of a bright compact radio source at 
this location. To obtain a better match with the observed ra- 
dio map, we reflected the map along the y-axis, i.e., placing 
the off-axis merger point on the other side (this is equivalent 
of changing the impact parameter from b to —6). The simu- 
lated profiles all peak around x = kpc, while the observed 
profile does not. From these profiles it is therefore difficult 
to constrain the impact parameter. However, by looking at 
Fig. [4] we find that b < 3r c ,i is required to obtain a reason- 
able match. For larger impact parameters the curvature of 
top shock does not match the observed one (i.e., this can be 
seen on the right side of the top relic in the R21b4 panel 
from Fig. 

4-1.3 Viewing angle 

We varied the viewing angle (i) by rotating our computa- 
tional box around the x-axis before projecting it along the 
line of sight. The extent of the emission towards the cluster 
center increases by about a factor of three from i — 5° to 
i = 30°, see Fig. [4] Radial profiles of the brightness distribu- 
tion are shown in Fig. [5] (middle panel). The outer edge of 
the northern relic in CIZA J2242. 8+5301 is very pronounced 
and the brightness of the radio emission drops also quickly 
towards the cluster center, therefore the relic is likely seen 
close to edge on, roughly under an angle < lOdeg. In all 
cases, the simulated brightness profiles do not reproduce the 
very quick drop towards the cluster center as it is observed 
for the northern relic. The southern relic is more diffuse and 
extended in CIZA J2242. 8+5301, as is the case in our sim- 
ulated maps. 



4-1-4 /3-model 

One of the parameters which could influence the size of the 
radio relics is the core radius used in our /3-profile (Eq.[7|. 
Because the total mass is fixed in our simulations the core 
radius is set by the j3 parameter in Eq.[7] We changed /3 from 
the default value of 2/3 to 0.5 and 0.75 for the most massive 
cluster to investigate the effect this has on the produced 
radio maps. For /? = 0.5, the core radius decreases from 
134 kpc to about 83 kpc, while for /3 = 0.75 the core radius 
increases to 163 kpc. We run the simulations with these 
different core radii and /3 values for the standard R21b0 
merger scenario. The resulting radio maps are displayed in 
Fig.© 

The change in r c ,i has a direct influence on the size 
of the resulting shock waves and radio relics formed. For 
r c ,i — 83 kpc, the northern and southern relics are almost 
a factor of 2 smaller in size. For r c ,i = 163 kpc, the relics 
have a slightly larger size and are similar to the R1.51b0 
merger (see Fig. [4|. This implies that our best fitting mass 
ratio found in Sect. |4. lTT] depend on the adopted core radii. 
In particular if the most massive cluster has a larger core 
radius this mimics a merger with a smaller mass ratio (e.g., 
1.5 : 1), see Fig. Bj By increasing the core radius of the 
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temperature -t density (z = 0) 300.09 Myr _ temperdture + density (z = 0) 900.08 Myr tempernture + density (z = 0) 1200.65 Myr 




Figure 1. Temperature and density evolution for a 2:1 merger with zero impact parameter. The color image displays the tempera- 
ture at a slice z = 0, while the white contours follow the density. The density contours are drawn at levels of log 10 (p([g cm -3 ]) = 
[-27.9, -27.8,-27.7,...]. 




Figure 2. Temperature and density evolution for a 2:1 merger with an impact parameter b = 4r c i. The color image displays the 
temperature at a slice z = 0, while the white contours follow the density. The density contours are drawn at levels of log 10 (p([g cm -3 ]) = 
[-27.9,-27.8,-27.7,...]. 
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Figure 3. Velocity field for a 2:1 merger with impact parameters 6 = kpc (left) and b = 4r Cj i (right) for the slice 2 = 0. The two 
snapshots were chosen such that the center of the two shock waves have a linear separation of 2.8 Mpc. The color image displays the 
temperature distribution. The maximum absolute velocity is ~ 1400 km s — 1 . 



Table 1. List of simulations with their defining parameters 



run 


mass ratio" 


impact parameter (b) 


r c ,l, r c ,2 


Tavg.n r avg,2 fy c comment 




Mi : M 2 


kpc 


kpc, kpc 


keV, keV km s" 1 



RllbO 


1:1 







115, 115 


4.1, 


4.1 


-937 




R1.51b0 


1.5:1 







126, 101 


4.5, 


3.7 


-927 




R21b0 


1:1 







134, 90 


4.7, 


3.3 


-922 




R2.51b0 


2.5:1 







139, 84 


4.8, 


3.1 


-926 




R31b0 


3:1 







145, 80 


5.0, 


3.0 


-944 




R21b0.5 


2:1 


lr c ,l 


= 67 kpc 


134, 90 


4.7, 


3.3 


-922 




R21bl 


2:1 


lrc,i 


= 134 kpc 


134, 90 


4.7, 


3.3 


-921 




R21b2 


2:1 


2r c ,i 


= 268 kpc 


134, 90 


4.7, 


3.3 


-916 




R21b3 


2:1 


3r c ,i 


= 402 kpc 


134, 90 


4.7, 


3.3 


-908 




R21b4 


2:1 


4r c ,i 


= 536 kpc 


134, 90 


4.7, 


3.3 


-897 




R21b5 


2:1 


5r c ,i 


= 670 kpc 


134, 90 


4.7, 


3.3 


-883 




R21b0/30.5 


2:1 







82.6, 90 


4.7, 


3.3 


-922 


/3i = 0.5 


R21b0/30.75 


2:1 







162.5, 90 


4.7, 


3.3 


-922 


ft = 0.75 


R21b0cc 


2:1 







134, 90 


4.7, 


3.3 


-922 


cool core'' for cluster 1 


R21b0ss (A=0.2, A = 100 kpc) 


2:1 







134, 90 


4.7, 


3.3 


-922 


substructure 


R21b0ss (A=0.3, A = 75 kpc) 


2:1 







134, 90 


4.7, 


3.3 


-922 


substructure 


R21b0ss (A=0.3, A = 150 kpc) 


2:1 







134, 90 


4.7, 


3.3 


-922 


substructure 


R21b0ss (A=0.3, A = 200 kpc) 


2:1 







134, 90 


4.7, 


3.3 


-922 


substructure 


R21b0ss (0.4, 200 kpc) 


2:1 







134, 90 


4.7, 


3.3 


-922 


substructure 



subscript 1 refers to the most massive subcluster, Mi + M2 = 5.5 X 10 14 
b T min ,i/T avgi i = 0.3, r cool ,i = 50 kpc and , a cool ,i = 2 
c starting velocity in y-direction for the least massive cluster 2 (Eq|l3) 
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R1.51b0 R21bO R2.51b0 R31bO 




R21b0, indinotion-5 deg R21bO, inclinolion— 1 deg R21b0, inclfnotion = 20 deg R21b0, inelinotion-50 deg 




Figure 4. Simulated radio emission for cluster mergers with different mass ratios (top row), impact parameters (middle row) and viewing 
angles (bottom row). The radio maps to the hydrodynamical snapshot where the two shock are separated by 2.8 Mpc. The images span 
3.4 X 3.4 Mpc. Blue contours show the observed emission at 1382 MHz from CIZA J2242. 8+5301 with the WSRT. The orange colored 
image and black contours display the synthetic radio image. Contours are spaced by factors of two in brightness 



less massive subcluster the size of the smaller southern relic 
can be increased, to some extent mimicking a 2.5 : 1 merger. 
Taking different core radii, the range of matching mass ratios 
for CIZA J2242. 8+5301 is therefore larger and lies roughly 
between 1.5 : 1 and 2.5 : 1 . 



4-1.5 Cool core 

In massive relaxed galaxy clusters the ICM temperature de- 
creases strongly towards the center. To investigate the ef- 
fect of such a cool core, we add a central cool core to the 
most massive subcluster using Eq. |10| We take typical value 
of Tmin/ravg = 0-3, r coo i = 50 kpc and r coo i = 2 (e.g., 
|Vikhlinin eFaLpOOo) and use the standard R21b0 merger 
scenario. We find that the addition of a central cool core 
does not influence the size and shape of the resulting radio 
relics. The resulting radio map looks almost identical to the 
default R21b0 merger map, see Fig. [4] 



4.2 Spectral index 

We computed a spectral index map between 610 and 
1382 MHz for the R21b0 merger, again taking a constant 
magnetic strength of 5 /iGauss, see Fig. [7] We find that the 
spectral index in front of the two relics is ~ —1-1, while 
the spectral index steepens quickly to ~ —2 in the direction 
of the cluster center. After the initial spectral steepening 
to a ~ —2, the spectral index remains constant. The radio 
emission in this region comes from (i) aged radio plasma 
with a steep spectrum in the post-shock region and (ii) pro- 
jected radio emission from the front of the shock wave with 
a flatter spectral index of ~ — 1. 

The observed spectral index in front of the northern 
CIZA J2242. 8+5301 relic is flatter, with a ~ -0.6 or -0.7, 
than in our simulated spectral index map. The observed 
spectral index corresponds to a M. = 4.6 shock. In the our 
R21b0 merger, we find a Mach number of about 2.7-2.4 
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R21b0ff0.5O R21 b0ff0.75 R21bOstop 




Figure 6. Top row: Simulated radio emission for R21b0 cluster merger with different values of f}\ and r G ,i and for a R21b0 cluster merger 
with the two halos sticking at core passage, which represents the extreme case of infinite tidal friction. Bottom row: Simulated radio 
emission for a R21b0 cluster mergers with substructure. Contours and colors are the same as in Fig. [4] The images span 3.4 X 3.4 Mpc. 



for the upper shock wave. This translates into a 1YI 



-0.8 



to —0.9. However, in the simulated map a is steeper by 
about 0.3 units at the front of the relic. This difference is 
caused by the limited spatial resolution (~ 50 kpc) of our 
synthetic radio map: radio emission directly at the front of 
the shock with aj n j w —0.85 is always mixed with some 
emission from particles that have already have undergone 
some synchrotron losses. For CIZA J2242. 8+5301 this im- 
plies that the true injection spectral index is probably even 
flatter than —0.6, if the shock front has a similar shape as 
in our simulations. 



5 SUBSTRUCTURE & CLUMPING 

Real mergers are more complex than in our idealized simu- 
lations. Hence, the detailed properties, such as the shapes or 
brightness profiles of radio relics, do not necessarily match 
with the observations. In particular the presence of less mas- 
sive substructures and galaxy filaments extending from the 
large scale cosmic web will modify the shock morpholo- 
gies and locations (e.g., Paul et al. 20101. Indeed, 



there 

are additional smaller relics located in CIZA J2242. 8+5301 
that points towards a more complex merger event. These 



more complex relics are also seen in simulations which in- 
clude large-scale structure formation with radio emission 



from shocks (e.g., Pfrommer et al.||2008 Hoeft et al. 2008 



Battaglia et al.|2009 \ . Although we cannot easily embed our 
simulation into a large-scale cosmological environment, we 
investigate the effect of small scale (< 200 kpc) substruc- 
tures/clumps. Matter falling onto clusters is expected to be 
clumpy in nature and the clumps have typical dimensions 



of galaxies or small galaxy-groups (Roncarelli et al. 2006) 



If the gas in the cluster's outskirts is clumpy, and this is 
not taken into account, it can lead to an overestimation of 
the baryon fraction and affect the derived entropy and pres- 
sure profiles (e.g., Simionescu et al.|2011 |. The presence of 
well-defined radio relics with sharp outer rims can be used 
to constrain the amount of clumping in the ICM at large 
distances (> 1 Mpc) from the cluster center. 

We modified our initial density and temperature pro- 
files by adding sinusoidal density fluctuations throughout 
the computational volume. The two parameters which con- 
trol this are the relative amplitude (A) and wavelength 
(A) of these density fluctuations. These sinusoidal waves 
were inserted parallel to the x, y and z-axes. To keep the 
pressure unchanged, we also included temperature fluctua- 
tions that are inversely proportional to the density fluctu- 
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Figure 5. Top: integrated brightness profiles across the y-axis 
for the top relic. Black dashed line displays the observed profile 
for the northern relic in CIZA J2242.8+5301.Solid lines, black: 
6 = kpc, purple: b = lr c , blue: b = 2r c , green: b = 3r c , orange: 
b = 4r c , red: b = 5r c . Middle: integrated brightness profiles across 
the x-axis from \x\ < 100 kpc for the top relic for relics seen under 
different viewing angles. Black dashed line displays the observed 
profile for the northern relic in CIZA J2242. 8+5301. Solid lines, 
black: i = 0°, purple: i = 5°, blue: i = 10°, green: i = 20°, red: 
i = 30°. Bottom: integrated brightness profiles across the x-axis 
from \x\ < 100 kpc for simulations including substructure. Black 
dashed line displays the observed profile for the northern relic in 
CIZA J2242. 8+5301. Solid lines, black: no substructure, purple: 
A = 0.2, A = 100 kpc, green: A = 0.3, A = 150 kpc, red: A = 0.4, 
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Figure 7. Simulated spectral index map between 610 and 
1382 MHz for the R21b0 merger. The small scale spectral in- 
dex variations seen in the post-shock region are due to the finite 
number of particles used to compute the radio map. 



ations, see Fig. [§] for an example. We carried out several 
runs changing the amplitude and wavelength of the fluctu- 
ations. Radio images from runs with A = 0.2,0.3,0.4 and 
A = 100, 150, 200 kpc, respectively, are shown in Fig. [6] 
Brightness profiles along the y-axis are displayed in Fig.[5] 
(bottom panel). 

The resulting radio relics are similar in shape to the 
canonical R21b0 merger map, but the relics are wider and 
display brightness fluctuations in the radio emission. The 
effect of this substructure is that the outer boundary of 
the relics is less well defined. For relics observed close to 
edge-on, projection effects average out some of the vari- 
ations in the strength of the radio emission. The north- 
ern relic in CIZA J2242. 8+5301 has a smooth appearance 
and well-defined outer rim. This probably implies that few 
clumps are present this location in the ICM. We find that 
the radio maps with A > 0.3 do not provide a good 
match to the observed northern relic in CIZA J2242. 8+5301. 
There is little difference between runs with A — 0.3 and 
A = 75, 150, 200 kpc. We can therefore put upper limits on 
the amplitude of density fluctuations of roughly 30%, for 
substructures with scales of < 200 kpc, at the location of 
the northern relic. This means that clumping appears to be 
much less than that inferred from Suzaku observations of 
Perseus with y/< n% > / < n e > 2 ~ 3 — 4 near the virial 
radius (e.g., Simionescu et al. 20111. The reason for this 



could be that the shock waves in CIZA J2242. 8+5301 are 
expanding into a region sufficiently far from filaments such 
that it is largely free from infalling substructure. 

A more realistic representation of density clumps could 
probably be used to put tighter constraints on the amount 
of substructure. This is beyond the scope of this paper and 
we leave this for future work. 
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Figure 8. Left: Temperature distribution at the slice 2 = for 
the R21b0ss (A=0.4, A = 200 kpc) merger. 



6 DISCUSSION 

6.1 Effect of dark matter dynamics 

Clearly, a gross simplification of our setup is the treatment 
of the gravity in the form of rigid potentials that move on set 
trajectories without interacting. In reality, the dark matter 
particles of the progenitor halos will interact in the course 
of the merger to form a new halo. In the course of this, 
the potential changes mainly due to two effects: i) dynamic 
friction and ii) tidal shocking. The first one is the most im- 
portant effect. After core passage, it will cause the halos to 
slow with respect to each other until they turn around on 
their orbits. As the merger shock waves are running out of 
the cluster faster than the speeds of the halos with respect to 
each other, this will not change the results dramatically. To 
test the magnitude of this effect, we ran a simulation where 
we made the two halos stick at core passage, which repre- 
sents the extreme case of infinite tidal friction. The result 
is shown in Fig. [6] The relics are now somewhat narrower 
(about 10% for the top relic) but apart from that there are 
not many changes in the radio map even though the tem- 
perature and density are quite different in the center. The 
main change are the lower Mach numbers (which drop from 
2.7 to 2.2, and from 2.2 to 1.9 for the top and bottom relics 
compared to the default R21b0 merger), but since this af- 
fects both for the top and bottom shocks the end result is 
not too different from the default R21b0 map. Finally, we 
ignore the effect of tidal shocking, which leads to heating of 
the central DM orbits of the cluster and leads to a puffing up 
of the potential. This will affect the cluster potential within 
the core radius but will leave the regions where the shocks 
are located largely unaffected. For a proper simulation of 
this effect, one would need to resort to a simulation with a 
dynamic DM component. 



6.2 Relic width and brightness profiles 

One of the most intriguing properties of the northern relic in 
CIZA J2242. 8+5301 is its very narrow width and rather uni- 
form luminosity along its extent. From our simulated radio 
map we find the brightest part of the top relic to be ~ 20% 
wider than the observed width. In Sect. 16.11 we found a 
smaller width by including infinite tidal friction. This affects 
the shock downstream velocity, which determines how far 
the the synchrotron emitting electrons are carried away from 
the shock front (e.g., 



Markevitch et al. 20051. The relic's 



width is linearly proportional to the downstream shock ve- 
locity. In case of infinite tidal friction the simulated radio 
relic is only ~ 10% wider compared to the observed width. 
In addition, the width depends on the adopted magnetic 
field strength, but increasing B to 7 /iGauss decreases the 
width only by a few percent. This implies that projection 
effects and the finite resolution of our simulations play and 
important role in determining the width. 

For the southern relic, the emission fades slowly towards 
the east and west ends. The reason why this does not hap- 
pen for the top relic in the simulation is because of the quick 
drop in the Mach number at \x\ > 850 kpc, see Fig. [3] (left 
panel). However, for the CIZA J2242. 8+5301 northern relic 
the integrated brightness profile across the y-axis (Fig[5] top 
panel) drops much quicker and is flatter than in the simula- 
tions. The reason for this is unclear, it could be caused by 
variations in the magnetic field strength, Mach number, or 
because the shock surface is more complex. The Mach num- 
ber found from the injection spectral index is higher than 
that in the simulations (4.6 versus 2.7 for the northern relic), 
this difference is puzzling. It could indicate that our temper- 
ature profiles are not correct or that Eq. [3j which links a.- m j 
and M, is not directly applicable. X-ray observations are 
crucial here to directly determine the Mach number from 
the temperature jump. 



6.3 Origin of single radio relics 

Another interesting outcome of our simulations concerns the 
existence of single radio relics, e.g., the relic in A521 (Gi- 



acintucci et al.|2008 l. In binary cluster mergers, two shocks 
should be generated and hence one expects the majority of 
relics to come in pairs. In the simulations with mass ratios 
> 3 very little radio emission is generated by the shock wave 
in front of the less massive subcluster. The Mach number of 
these shocks is quite low (A4 < 2) and the shock surfaces are 
small, which results in very little observable radio emission. 
Based on the simulations, one would thus expect most dou- 
ble relics to be located in merging systems with small mass 
ratios, while single should be found in mergers with larger 
mass ratios. 



6.4 A quantitative metric for the goodness of fit 

We did not define a quantitative metric for the goodness 
of fit between the simulated and observed radio images. 
This paper is meant to be illustrative of the potential of 
the method and in this work we chose to apply it only to 
the relics in CIZA J2242. 8+5301. We also note that the po- 
sition, shape, and surface brightness distribution of relics 
are influenced by the cluster environment (e.g., connected 
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galaxy filaments) which affect the derived merger param- 
eters. However, for future work when constraining merger 
parameters for larger numbers of double relics systems, it is 
important to develop a metric for the goodness of fit. We 
discuss here a possible approach and address some of the 
issues involved. 

A first step should be the removal (or masking) of point 
sources and extended radio galaxies. In addition, one may 
have to subtract the emission from a radio halo in the cluster 
center. Removing extended radio galaxies and a radio halo 
is not entirely trivial as the spatial scales of these structures 
can be comparable to that of relics. However, this should 
still be feasible as double radio relics are relatively bright 
and their morphology is quite different from radio halos or 
radio galaxies. 

For constraining the merger parameters there are two 
aspects that can be considered: (i) the spatial information 
of the location of the relics and (ii) the brightness distri- 
bution. We ignore information about the spectral and po- 
larization properties because this complicates matters and 
this information is not available for all double relics. Due to 
the uncertainties in the magnetic field strength the bright- 
ness of the emission is probably not a very good criteria 
for a quantitative metric. It is also not entirely clear under 
what conditions relics precisely form, i.e, in the case of Abell 
2146 two shocks are found in a Chandra X-ray image, but no 
bright radio relics are detected ( |Russell et al.|2010|[20Tl) . 

Radio relics, especially their outer edges, should directly 
trace the location of the merger shock waves and thus pro- 
vide most of the constraints for the merger scenario. This 
is also what we mainly used in our qualitative comparison 
for the relics in CIZA J2242.8+5301. The viewing geometry 
is mostly constrained from the downstream extent of the 
relics. In the light of what we discussed above, we propose 
a two-dimensional spatial cross-correlation. This should be 
repeated for every possible rotation angle for a given image 
from the simulations. The height of the peak in the cross- 
correlation then gives a quantitative metric for the goodness 
of fit for a single simulated image. Image scaling should not 
be considered here as this relates to the absolute sizes and 
distances of the relics which evolve over time, and the time 
since core passage is one of the parameters we wish to con- 
strain. The process can then be repeated for different merger 
parameters (mass ratios and impact parameters), snapshots 
in time and viewing angles to find the best fit. Because the 
brightness distribution is uncertain, we propose to apply a 
fixed cut in the radio flux (a few times the noise level in the 
observed radio image), and put all pixels values to either 1 
(emission) or (no emission) . This is done also for the simu- 
lated images. The overall normalization of the simulated ra- 
dio emission can be fixed by normalizing to the total flux in 
the observed image. Instead of a cross-correlation, one could 
also consider a simple image subtraction (after x, y-shifting 
and rotation) followed by taking the sum of the absolute dif- 
ference. It may be worth giving the outer edges of the relics 
a somewhat large weight in the process, for example by set- 
ting these to a higher value (e.g., 2), some experimentation 
will be needed here. Applications of these cross-correlations 
shall be reserved for future work. 



7 SUMMARY 

We have presented simulated radio maps of idealized binary 
cluster mergers, varying the mass ratios, impact parameters, 
viewing angles and density/temperature profiles. The radio 
emission in these simulations is produced by relativistic elec- 
trons which are accelerated via the DSA mechanism in the 
shock waves created by the merger event. 

We have used these simulated radio maps to constrain 
the merger scenario for the cluster CIZA J2242. 8+5301, 
which hosts a Mpc scale double radio relic. We find that a 
merger with a mass ratio between 1.5 : 1 and 2.5 : 1 provides 
the best match to the observed radio emission in this cluster. 
The main uncertainty in the derived mass ratio comes from 
the unknown core radii of the clusters before the merger. 
The viewing angle is mostly constrained by the extent of 
the top relic towards the cluster center, additional substruc- 
ture/clumps also increase the relic's extent in the down- 
stream region. Since the observed extent is only 55 kpc this 
implies that the relic is seen close to edge-on (i.e., i < 10°). 
The impact parameter is difficult to constrain as the ob- 
served brightness profile along the x-axis is not reproduced 
for any our merger scenarios. Only for b > 3r c ,i « 400 kpc 
we find that the curvature of the simulated relics do not 
match with the observations. From our simulations we find 
that the time since core passage does not depend much on 
the adopted initial conditions. We therefore constrain this 
number to 1.0 Gyr, although it also depends on the adopted 
total mass of the system and temperature profile. 

One of the main uncertainties for comparing the sim- 
ulated radio emission with observations, in particular the 
brightness profiles, is the structure and strength of the mag- 
netic field. Especially the brightness profile of the main relic 
along the relic's length is not well reproduced by our sim- 
ulations. Little is known about the magnetic field strength 
at distances larger than 1 Mpc from the cluster center. Us- 
ing equipartition arguments magnetic fields strengths of the 
order of ~ 1 /iG have been derived (e.g 



Govoni & Feretti 



2004|. Upper limits on the IC emission from radio relics in- 



dicate magnetic fields of a few /xG (e.g., Nakazawa et al. 
|2009| |Finoguenov et al. 20101. Faraday rotation measures 
give similar results (e.g., | Bonafede et al.|[2010[ ). However, 
the spatial variations of the magnetic field strength at these 
distances from the cluster center are unknown. In addition, 
several smaller relics are found in the cluster which are not 
present in our simulations. These relics are probably caused 
by additional substructures which are not included in our 
initial conditions. 

We find that radio relics can be used to put constraints 
of the amount of clumping or substructure in the ICM at 
large distances (~ 1 - 2 Mpc) from the cluster center, under 
the assumption that relics are tracers of shocks. The sharp 
outer boundary of the northern relic in CIZA J2242. 8+5301 
implies density fluctuations with an amplitude < 30% of the 
average density (the size of fluctuations being ~ 200 kpc). 

The simulations also indicate that most double relics 
should be located in merging system with mass ratios be- 
tween 1 : 1 and 1:3, while single radio relics from in clusters 
that undergo less energetic mergers with larger mass ratios. 

We conclude that double radio relics can be used to 
constrain cluster merger events, which is especially useful 
when detailed X-ray observations are not available. Relics 
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are often located at large distances from the cluster center, 
where the density is so low that even deep X-ray observa- 
tions cannot be used to probe the physical conditions there 
(e.g., to measure the temperature jump). In addition, relics 
can be located in lower-mass systems where the thermal X- 
ray emission is not detected in the ROSAT All Sky Survey, 



for example ( Voges et al.||1999 1. Some examples of possible 
radio relics in poor systems without a X-ray counterparts 



have been reported by Brown & Rudnick ( 2009 1; van Weeren 
eFHTl p009l [20TT I . 



One important parameter is the strength and structure 
of the magnetic field at the location of the radio relics. MHD 
simulations of cluster mergers will be needed to study the 
variations and strength of the magnetic field at large dis- 
tances from the cluster center, and investigate the effect this 
has on the derived merger parameters from radio maps. 
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